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ABSTRACT 



Aims. We search for galaxy counterparts to damped Lyman-ff absorbers (DLAs) at z > 2 towards nine quasars, which have 14 DLAs 
and 8 sub-DLAs in their spectra. 

Methods. We use integral field spectroscopy to search for Lyo- emission line objects at the redshifts of the absorption systems. 
Results. Besides recovering two previously confirmed objects, we find six statistically significant candidate Lya emission line objects. 
The candidates are identified as having wavelengths close to the DLA line where the background quasar emission is absorbed. In 
comparison with the six currently known hya emitting DLA galaxies the candidates have similar line fluxes and line widths, while 
velocity offsets between the emission lines and systemic DLA redshifts are larger. The impact parameters are larger than 10 kpc, and 
lower column density systems are found at larger impact parameters. 

Conclusions. Assuming that a single gas cloud extends from the QSO line of sight to the location of the candidate emission line, we 
find that the average candidate DLA galaxy is surrounded by neutral gas with an exponential scale length of ~5 kpc. 

Key words, galaxies: formation - galaxies: high redshift - galaxies: quasars: absorption lines 



1. Introduction 

Galaxy counterparts to Damped Lyman-a systems (DLAs) seen 
in quasar (QSO) spectra have been suggested to be (proto)-disk 
galaxies with line of sight clouds of neutral gas with column 
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densities JV(Hi) > 2 x 10 20 cnr 2 dWolfe et al.l ll986). Analyses 
of absorption line profiles have indicated that rotational compo- 
nents with velocities of ~200 km ~' can be involved in these 
systems suggesting that DLAs reside in large disk galaxies 
dProchaska & Wolfel \l99% iLedoux et ai1ll998al) . On the other 
hand, numerical simulations show that in a hierarchical forma- 
tion scenario merging proto- galactic clumps can also give rise to 
the observed line profiles dHaehnelt et alj[l998h . 

A large fraction of the neutral hydrogen pres ent at z > 2 is 
conta i ned in high column density DLA systems ( Lanzetta et alj 
Il995t IStorrie-Lombardi & Wolfel l2000t iPeroux et alJl2001h . In 
addition to the classical DLAs, clouds with column densities 
10 19 < N(Ui) < 2 x 10 20 cm~ 2 also show some de- 
gree of damping wings, which is characteristic of DLA sys- 
tems. It is suggested that sub-DLA systems cont ain a signifi- 
cant f raction of the neutral matter in the Universe dPeroux et alj 
12003b . Metallicity studies have shown that the properties of 
the sub-DLA systems are simil ar to those of DLA systems 
dDessauges-Zavadskv et alj|2003l) . ap art from the latter category 
havin g large ionisation corrections dProchaska & Herbert-Fortl 
2004). 



Send offprint requests to: L. Christensen 

* Based on observations collected at the Centra Astronomico 
Hispano Aleman (CAHA), operated by the Max-Planck Institut fiir 
Astronomie and the Instituto Astrofisica de Andalucia (CSIC). 



The association of DLAs with galaxies has been a subject of 
much study. Originally, either space-based or ground-based deep 
images we re obtained to identify objects near the l i ne of sight to 
the Q SOs dSteidel et alii 19951. iLe Brun et al]ll997tlWarren et alJ 
2001). To confirm nearby objects as galaxies that are respon- 
sible for the DLA lines in the QSO spectra, follow-up spectra 
are needed to find the galaxy redshifts. At z < 1, confirmations 
of 14 systems exist to date {Rao et alj 120031: IChen & Lanzetta! 
l2003l:lLacv et al.l2003tlChen et al.l2005l and references therein), 
while at z > 2 only 6 DLA galaxies are confirmed through spec- 
tros copic observations of Lyq emission from the DLA galax- 
ies dM0ller & Warrenll993l:lDiorgovski et alJl996tlM0ller et all 
U998l:lLeibundgut& Robertsonll 19991: lM0ller et alj|2002l 12004). 
three of which are located at the same redshifts as the QSOs 
themselves. Other technique s to identify DLA galax ies involve 
narrow-band imaging (e.g. Fvnb o et ail 1 19991 l2000h or Fabry- 
Perot imaging. A Fabry-Perot imaging study of several QSO 
fields did not result in de tections of emission from DLA galax- 
ies dLowenthal et al.|[l995l) . while recently the sa me method was 
used to identify a few emission line candidates (Kulkar ni et al.l 
120061) . 

Integral field spectroscopy (IFS) presents an alternative that 
provides images and spectra at each point on the sky simultane- 
ously. This technique can be used to look for emission line ob- 
jects at known wavelengths, but unknown spatial location. This 
technique is ideally suited to look for Lya emission lines from 
the galaxies responsible for strong QSO absorption lines. At the 
Lya wavelength corresponding to the redshift of the DLA sys- 
tem, the QSO emission has been absorbed, enabling us to lo- 
cate emission line objects very near to the QSO line of sight. 
Because of the large column densities in DLAs and the res- 



2 



L. Christensen et al.: An IFS survey for high-z DLA galaxies 



onant nature of Lya photons the corresponding emission line 
may be offset in velocity spac e relative to the DLA line (e.g. 
Leibundgut & Robertson 1999), but such an offset is not always 
observed (e.g. |M0ller et alj|2004l) . 

IFS searches for emission from DLA galaxies towards 
two QSOs have re s ulted in upper limi ts for their fluxes 
dPetitiean et al.l Il996t iLedoux et al]|1998bl) . while a sub-DLA 
galaxy pr eviously known to be a Lya emitter was confirmed 
with IFS (Christens en et al.l l2.004t). Here we present a survey us- 
ing IFS to look for Lya emitting DLA galaxies. Section [2] de- 
scribes the sample of QSOs included in the survey, which are 
known previously to have DLAs and sub-DLAs in their spec- 
tra. Section [3] describes the observations and data reduction. In 
Section [4] the method of detecting emission line candidates is 
described. Section[3]presents the results and comments on each 
object. Properties of the Lya emission candidates detected in 
the survey in relation to the six previously known Lya emitting 
DLA galaxies are presented in Section [6] Section [7] summarises 
our findings. A flat cosmology with Ho - 70 km s _1 Mpc -1 , 
Q m = 0.3, and £2a = 0.7 is used throughout. 

This study, as well as previous ones that try to identify the 
host galaxies of DLA systems, can be biased since the galaxy 
observed at the right redshift likely belongs to the brightest emis- 
sion line object close to the line of sight. In the case that the 
host galaxy is a much fainter galaxy in a group, it will not be 
identified correctly. In the remaining part of the paper, an 'iden- 
tified' DLA galaxy refers to observations that show (line) emis- 
sion from independent observations, while the 'candidates' are 
only reported in these IFS observations. Although extensive tests 
are done on the data to distinguish the candidates from potential 
artifacts, independent observations are needed to prove them as 
Lya emitters connected with the DLAs. 

2. Sample selection 

We selected a number of DLA systems without previous de- 
tections of associated Lya emission. The selected QSOs with 
known DLAs were chosen based on the following criteria 

1. JV(Hi)>2x 10 20 cnr 2 

2. DLA redshift (2 < z < 4) 

3. Northern hemisphere object 

The first criterion includes only classic DLA systems. Many 
QSO spectra show additional sub-DLA systems. Although their 
relationship to DLAs is still debated, these systems were in- 
cluded in the survey because of their probable physical associa- 
tion with galaxies. 

To increase the sample size with a minimum number of 
pointings we preferentially selected QSOs with multiple DLAs. 
IFS covers a range of wavelengths, and correspondingly Lya 
emission at a large range of redshifts in the line of sight for each 
QSO. However, in retrospect, this can affect the emission line 
detections, because extinction in foreground DLAs could ob- 
scure emission from background ones when the galaxies lie in 
the same line of sight. Hence, upper limits on detections of the 
higher redshift systems can be biased. 

From the list of DLA systems compiled by S. CurrarQ, 
we found 66 QSOs matching these criteria in 2003. More re- 
cently, detections of DLAs in the Sloan Digitized Sky Survey 
QSO spectra have greatly increased the number of know n DLAs 
dProchaska & Herbert-Forl2004t iProchaska et al.1120051) . A sys- 
tematic survey of all 66 objects would require a large amount 

1 http : //www . phys . unsw . edu . au/~s j c/dla/ 



of time with present instruments, so we selected a few systems 
based on their observability during the allocated observing runs. 
We avoided DLAs with Lya absorption lines close to sky emis- 
sion lines. 

The total sample consists of 9 QSOs with a total number of 
14 DLA systems as listed in Table Q] These QSOs have an ad- 
ditional 8 sub-DLAs which are included in the survey. Because 
of the small number of DLAs involved in the survey, a proper 
statistical study is not the aim of this paper. Instead we focus 
on a few systems to exploit the benefits of IFS for this kind of 
investigation. 

To study the applicability of IFS in identifying DLA galax- 
ies we initially observed DLA galaxies where Lya emis- 
sion had been reported previously in the literature. Two of 
these systems could be observed during o ur runs; Q2233+13 1 
and PHL 1222 origin ally i dentified by ISteidel et all d 1995b . 
iDiorgovski et alj d 19961) and lM0ller et all dl998l). Both objects 
are reported to have ext ended Lya emission dFvnbo et alJ fl999: 
IChristensen et al.l l2004). Table Q] includes these two previously 
known Lya emitting DLA and sub-DLA galaxies, although the 
criteria listed above are not satisfied. The absorption system to- 
wards Q2233+131 has a column density that classifies it as a 
sub-DLA. Unless otherwise noted, these two objects are kept 
separate from the detection of candidate emission line objects in 
the remainder of the paper. 

Most of the DLAs in the IFS study lie towards bright QSOs 
(R < 19). This ensured that the PSF variations as a function of 
wavelength could be determined, which was necessary for the 
subtraction of the QSO emission. Bright QSOs had larger resid- 
uals from the subtraction of the continuum emission, which po- 
tentially affected out ability to recover emission line objects that 
were offset in velocity space and located closer than 1" to the 
QSO line of sight. However, tests with artificial objects showed 
that this was a minor problem for the data set (see Sect. 15.31 ). 



3. Observations and data reduction 

Using the Potsdam Multi Aperture Spectrophotometer (PMAS) 
mounted on the 3.5m telescope at Calar Alto we observed the 
objects listed in Table [2] during several runs from 2002-2004. 
The PMAS integral field unit (IFU) was used in the standard 
configuration where 256 fibres are coupled to a 16x16 element 
lens array. During the observations each fibre covered 0'.'5x0'.'5 
on the sky giving a field of view of 8"x8". Each fibre resulted in 
a spatial element (spaxel) represented by a single spectrum. The 
256 spectra were recorded on a 2kx4k CCD which was read out 
in a 2x2 binned mode. With a separation of 7 pixels between in- 
dividual spectra, the fibre to fibre cross-talk was negligible (less 
than 0.4% for an extraction of all 7 pixels). Detailed overviews 
of the PMAS instrument and capabilities are given in Rot h et ail 
d2000ll2005h . 

For individual exposures a maximum time of 1 800s was used 
because of the large number of pixels affected by cosmic ray hits. 
Furthermore, because of varying conditions such as the atmo- 
spheric transmission and seeing, the total exposure time for each 
object was adjusted, or sometimes an observation was repeated 
under better conditions. The photometric conditions during ob- 
servations were monitored in real time with the PMAS acquisi- 
tion and guiding camera (A&G camera) which is equipped with 
a lkxlk CCD. Seeing values listed in Table [2] refer to the see- 
ing FWHM measured in the A&G camera images. Determining 
actual spectrophotometric conditions requires monitoring of the 
extinction coefficients which cannot be determined from the 
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Coordinate name 


Alt. name 


■-em 


Zabs 


1 net Affl-J T \ 

lOgiv^rl 1 ) 


Lre/hLJ 


rc; /TTl 


Kei. 


Q0151+048A 


PHL 1222 


1.93 


1.934 


20.36+0.10 






(i) 


Q0953+4749 


PC 0953+4749 


4.457 


3.404 
3.891 
4.244 


21.15+0.15 
21.20+0.10 
20.90+0.15 


>-2.178 
>-1.712 
-2.50+0.17 


>-2.09 

>-1.60 
-2.23+0.15 


(2,3) 


Q1347+112 




2.679 


2.471 
2.05 


20.3 
20.3+ 






(4,5,6) 
(7) 


Q 1425+606 


SBS 1425+606 


3.163 


2.827 


20.30+0.04 


-1.33+0.04 


>-1.03 


(8,9,10,21) 


Q1451 + 1223 


B1451 + 123 


3.246 


2.469 
3.171 
2.254 


20.39+0.10 
19.70+0.15 
20.30+0.15 


-2.54+0.12 
-1.87+0.16 
-1.47+0.17 


-1.95+0.16 
-1.62+0.15 
>-0.40 


(11,24) 
(19) 

(6,12,19) 


Q1759+7539 


GB2 1759+756 


3.05 


2.625 
2.91 


20.76+0.10 
19.8 


-1.21+0.10 
-1.65+0.01 


-0.82+0.10 
-1.26+0.01 


(13,18,20) 
(18) 


Q1802+5616 


PSS J1802+5616 


4.158 


3.391 
3.554 
3.762 
3.811 


20.30+0.10 
20.50+0.10 
20.55+0.15 
20.35+0.20 


-1.54+0.11 
-1.93+0.12 
-1.82+0.26 
-2.19+0.23 


>-1.55 
>-1.82 
>-1.74 
-2.04+0.22 


(23) 


Q2155+1358 


PSS J2155+1358 


4.256 


3.316 
3.142 
3.565 
4.212 


20.55+0.15 
19.94+0.10 
19.37+0.15 
19.61+0.15 


>-1.68 
-2.21+0.21 
<-2.40 
-2.18+0.25 


-1.26+0.17 
-1.85+0.12 
-1.27+0.16 
-1.92+0.11 


(3) 

(14,19) 
(14,19) 
(14,19) 


Q2233+131 




3.295 


3.153 
2.551 


20.0 
20.0 


>-l. 4+0.1 


>-1.04 


(15,16,17) 
(12,16) 



Table 1. List of the observed DLA and sub-DLA systems with column densities and metallicities taken from the litera- 
ture. * denotes a system where the reported jV(Hi) need s to be con firmed through high re solut ion spectroscopy . Refere nces 
for either DLA re dshifts or metalhc i ties: ( 1) iMejller et all (Il998l) (2) ISchneider et all dl99ll). ( 3) Prochask a et alj J2003d), (4) 
Smit h etal.1 (1 19861). ( 5) IWolfeetalJ dl986h (6) iTurnshek et al l dl989n. (7) IWolfeetalJ dl995h. (8) IChaffee et all d 19941). (9 ) 
Step anian et al.ldl996l). ( lO) lProchaska et al l d2002al) . (1 DlBechtoldld 19941) . (12 )lLanzetta et al.ld 19911). ( 13) lProchaska et al.l d20Ql), 
(14) |Peroux et all d2003l). (15)ISteidel et aiT d 19951). (16) iLu et all d 19971) . (1 7) EProchaska et al .1 d2003al). (18)lQutram et ail Jl 999). 
(19) iDessauges-Zavadsky etal.1 d2003l). (20) [Prochaska et al.l d2002bl). (21) iLu et al.1 dl996Th (22) iProchaska et al.1 (l2003cb . (23) 
Prochask aTt al.ld2003bl) . (24) lPetitiean et all (120001) . 



A&G camera images. In Table [2] 'stable' means that the pho- 
tometry of the guiding star was constant within 1% during the 
observations. 

The data were obtained using 2 gratings; one with 300 lines 
mrrT 1 and one with 600 lines mirT 1 , set at a chosen tilt to cover 
a selected wavelength range. The FWHM of the sky lines were 
measured to be 6.4 and 3.2 A, respectively. Observations of spec- 
trophotometric standard stars were carried out at the beginning 
and end of each night at the grating position used for the obser- 
vations. 

Data reduction was done by first subtracting an average bias 
frame. Before extracting the 256 spectra most c osmi c ray hits 
were removed by the routine described in lPvchl d2004l) . A high 
threshold was chosen such that not all cosmic rays were re- 
moved, because a low threshold would also remove bright sky 
emission lines from some spectra. Remaining cosmic rays were 
removed from the extracte d spectra using the program L.A. 
Cosmic dvan Do kkum 20011). 

The locations of the spectra on the CCD were found from 
exposures of a continuum lamp, taken either before or after the 
science exposures, using a tracing algorithm de veloped for th e 
IDL based PMAS data reduction package P3D dBeckerll2002l) . 
The spectral extraction was done in two ways; a 'simple extrac- 
tion' that added all flux from each spectrum on the CCD (i.e. an 
extraction width of 7 pixels), and another method that took into 
account the profile of the spectrum on the CCD. This second 
method assumed that the spectral profiles were represented by 
Gaussian functions (Gaussian extraction) where the widths were 



allowed to vary with wavelength. Widths were determined by fits 
to each of the 256 spectra as a function of the wavelength, and 
the extraction used these width in combination with the centre 
found from the tracing algorithm. The Gaussian profile is an ap- 
proximation because the profiles are not strictly Gaussian. The 
second method increased the signal-to-noise ratio by >10% for 
faint objects and therefore unless otherwise noted, the results 
from the 'Gaussian extraction' data cubes will be reported (see 
also ISanchezll2006l) . 

After extraction, the spectra were wavelength calibrated us- 
ing exposures of emission line lamps taken just before or after 
the observations. The wavelength calibration was done using the 
P3D reduction tool. Comparisons with sky emission lines indi- 
cated an accuracy of the wavelength calibration of about 10% of 
the spectral resolution. 

The spectra show a wavelength dependent fibre to fibre trans- 
mission. To correct for this effect, we extracted sky spectra ob- 
tained from twilight sky observations in the same way as the sci- 
ence observations. A one dimensional average sky spectrum was 
calculated. Each of the 256 spectra were divided by this average 
spectrum, and the fraction was fit by a polynomial function to re- 
duce noise. These polynomials were used to flat field the science 
spectra. 

Before combining individual frames, the extracted spectra 
were arranged into data cubes. Each data cube was corrected 
for extinction using an av erage extinction curve for Calar Alto 
dHopp & Fernandezl2002f) . The data cube combination took into 
account a correction for the differential atmospheric refraction 
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QSO 


date 


exp.time 


grating 


A coverage 


seeing 


conditions 






(s) 




(A) 










(^U1M+U48A 


1Af\o AO T7 

2003-08-27 


C , . 1 o A A 
DX LOW 


V600 


3500- 


-5080 


0.8- 


-1.2 


stable 


Q0953+4749 


ZUU4-U4- 1 o 


A \ t 1 O A A 

4x lqxjO 


V300 


3630- 


-6980 


0.9 




stable 




ZUU4-U4-Z 1 


1 onn 
JX loUU 


V300 


3630- 


-6980 


1.0 




non-phot. 


Q1347+1 12 


2004-04-20 


7x 1800 


V300 


3630- 


-6980 


0.6 




non-phot. 




1 AA/1 A A If! 

2UU4-U4-1V 


. . i O AA 


V300 


3630- 


-6750 


1.0 




stable 


Q1451 + 1223 


-1 A /I 1^7 

2004-04-17 


^7^,1 O f\f\ 

7x 1800 


V300 


3630- 


-6980 


0.8 




non-phot. 


nncn . neon 

Q1759+7539 


1AA /I A/1 11 


-t, , i onn 

/X 1800 


V300 


3630- 


-6980 


1.0- 


-1.5 


non-phot. 


1^1802+5616 


^nno a/; 1 O 

2003-06- Is 


2X 1800 


V600 


5100- 


-6650 


1.0 




non-phot. 




onm on 
2003-06-20 


3X toOO 


V600 


5100- 


-6650 


1.0 




non-phot. 




onm n/; o i 
2003-00-21 


a w 1 onn 
4X 1800 


^ 7 /inn 
V600 


5100- 


-6650 


1.8 




non-phot. 




2003-06-22 


6x1800 


V600 


5100- 


-6650 


0.9 




stable 


Q2155+1358 


2003-08-26 


7x1800 


V600 


4015- 


-5610 


0.7 




stable 




2003-08-27 


4x1800 


V600 


4015- 


-5610 


0.8 




non-phot. 


Q2233+131 


2002-09-03 f 


4x1800 


V300 


3930- 


-7250 


1.0- 


-1.3 






2003-08-24 


6x1800 


V600 


4000- 


-5600 


0.6 




stable 




2003-08-25 


4x1800 


V600 


4000- 


-5600 


0.7 




non-phot. 



Table 2. Log of the observations. The last two columns show the average seeing during the inte grations and the photom etric 
conditions derived from the A&G camera images. 1 Results from these observations are published in Chris tensen et alj d2004t) . 



using a theoretical prediction (Filippenko 1982). Relative spa- 
tial shifts between individual data cubes were determined from 
a two-dimensional Gaussian fit to the QSO PSF at a wavelength 
close to the strong DLA absorption lines. 

Subtraction of the sky background was an iterative process 
because the locations of the emission line objects of interest were 
not known beforehand. PMAS, in the configuration used, does 
not have specifically allocated sky fibres. Instead, we selected 
10-20 fibres uncontaminated by the QSO emission and the av- 
erage spectrum was subtracted from all 256 spectra. Different 
spaxel selections were examined visually to select an appropri- 
ate sky spectrum which had no emission line or noisy pixels in 
the spectral region of interest. 

Flux calibration was done in the standard way using obser- 
vations of spectrophotometric standard stars. A one-dimensional 
spectrum of the standard star was constructed by co-adding flux 
from all 256 spaxels. This was used to create a sensitivity func- 
tion that could be applied to each of the 256 spectra in the science 
exposures. For non-photometric nights the flux calibrated spec- 
tra were compared with QSO spectra from the literature to esti- 
mate photometric errors. However, no correction factor was ap- 
plied to our spectra, because an intrinsic variability of the QSOs 
would make such scaling uncertain. For some cases we note in 
Sect. |5.5| that there are differences which could be caused by ei- 
ther non-photometric conditions or intrinsic variability. 

For reference we show spectra of the target QSOs in Fig.Q] 
Where present, metal absorption lines corresponding to the 
highest column density DLAs are indicated. For QSOs with 
multiple DLAs lines only the DLA lines and their redshifts 
are indicated since the wavelength coverage does not include 
lines redwards of the QSO Lya line. A detailed analysis of 
metal absorption lines req uires higher resolution spectroscopy as 
presented elsewhere (e.g. [Prochaska et alj 12003d : iPeroux et alj 
120031: iDessauges-Zavadskv et al.ll2003l) . DLA redshifts derived 
from the metal absorption lines were consistent with those re- 
ported in the literature within the accuracy of the wavelength 
calibration of the data cubes. 



4. Search for DLA optical counterparts 

The observations covered the wavelengths of Lya for all but one 
of the strong absorption systems listed in TableQ] Only the high- 
est redshift sub-DLA system towards Q2155+1358 was not cov- 
ered, i.e. the total number of systems included in this analysis is 
21 DLA and sub-DLA systems. 

4.1. Expected sizes 

For this project we are only interested in small wavelength re- 
gions corresponding to Lyor at the DLA redshifts, and thus the 
search for candidate galaxies could be carried out using cus- 
tomised narrow-band filters. IFS, on the other hand, has the ad- 
vantage that the widths of the narrow-band filters can be ad- 
justed to match those of the emission lines. Typically customised 
narrow-band filters have a larger transmission FWHM than the 
spectral resolution of IFS data. Hence, IFS allows detection of 
emission line objects with a higher signal-to-noise ratio than 
that possible with narrow -band filters. A disadvantage is the rela- 
tively small field of view of current IFUs, but this is not a serious 
conc ern. One ca n estimate the expected sizes of DLA galaxies 
(see IWolfe et ai1ll986l) . Using a Schechter luminosity function 
and a power-law relation between the disk luminosity and gas 
radius given by the Holmberg relation R/R* = (L/L*~f, one can 
calculate the expected impact paramet er. Combining f3 = 0.26 
found for DLA galaxies at z < 1 dChen et al.ll2005l) wit h the 
luminosity function in for z = 3 galaxies dPoli et al.ll2003l) one 
finds R* * 30 kpc. If DLA galaxies are similar to or fainter than 
L* galaxies this implies that DLA galaxies at z > 2 are expected 
to lie closer than ~4" from the QSO line of sight. The small 
field of view of IFUs is therefore well suited to search for Lya 
emission from DLA galaxies. 

The estimated galaxy sizes are highly dependent on the 
parameters of the DLA galaxy luminosity and slope [}. Most 
probably, high redshift DLA galaxies are not regular disks like 
those in the local universe. Numerical models of DLAs predict 
that the galaxies are mostly smaller than 10 kpc, while obser- 
vations that give limits on the star-formation rates associated 
with DLAs suggest that DLAs are located in neutral gas around 
Lyman break galaxies (Wolfe & Chen 2006). As DLA galaxies 
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at z > 2 are generally f ound to be fainter than an L* galaxy 
dColbert & Malkar3l2002l) . we choose to consider only objects 
with impact parameters smaller than 30 kpc for a more detailed 
analysis. 

The impact parameters that we measure in the data corre- 
spond to the radially projected distances so the real distances to 
the absorber can be larger. Two candidates are found at impact 
parameters larger than 30 kpc, and they are likely not associated 
directly with the absorbers themselves. 

4.2. Candidate selection 

Some Lya emission lines from DL A galaxies are offse t from the 
QSO- DLA line by ~ 200 km s" 1 (iMoller et al.ll2002l) . whereas 
Lya emission from high redshift galaxies can have even large r 
offsets from the galaxy systemic redshift (Sha plev et al.l l2003). 
We therefore chose to focus on regions in the data cubes with 
velocities ranging from approximately -1000 to +1000 km s _1 
from the DLA lines. 

First, the reduced data cubes were stacked in a two- 
dimensional frame and inspected visually around the DLA lines 
for emission line objects. When the spatial offset from the QSO 
is larger than the seeing, or alternatively when the QSO is very 
faint, emission line objects can be identified directly because 
of the ordering of the spectra in the stacked spectrum. Where 
no objects could be detected visually further sampling of the 
data cubes was necessary to increase the signal-to-noise ratio 
to detect candidate emission line objects. Inspections o f the data 
cubes was done using the Euro3D visualisation tool dSanchezI 
l2004h . 

From the reduced, sky-subtracted and combined data cubes, 
narrow -band images were created with an initial width of 10- 
15 A depending on the spectral resolution of the observations. A 
set of images was created offset by -10 to +10 A from the DLA 
line to allow for possible velocity shifts of the Lya emission line, 
and inspected visually for objects brighter than the background. 
If detected, spectra from these brighter regions were co-added 
and inspected for emission lines at the wavelength chosen in the 
narrow-band image. This step was necessary to discriminate be- 
tween emission lines and individual noisy spectra. It is known 
that three blocks of 16 fibres, i.e. 48 fibres in an area of 1'.'5 to the 
west in the field of view, have lower than average transmission. 
The effect of correcting for the total throughput was that these 
spectra had lower signal-to-noise ratios. When narrow-band im- 
ages were created from the cubes, the higher variance in these 
spaxels could result in extreme values, seemingly inconsistent 
with the neighboring spaxels. Only by looking at the spectrum 
associated with a bright spaxel could it be determined if an emis- 
sion line was present, or if the spectrum was just noisy. If an 
emission line was seen, a second pass narrow-band image was 
created using the value of the emission line width to increase the 
signal of the detection. A second pass one-dimensional spectrum 
was created after inspecting the narrow-band image for more 
bright spaxels surrounding the emission line candidate. This pro- 
cedure was iterated until the signal in either narrow -band images 
or spectra did not increase. We found that an interactive visual 
identification of faint emission lines was more effective than an 
automatic routine. 

To allow a better visual detection of emission line objects, 
the narrow-band images were interpolated to pixel scales 072 
pixel 1 as shown in Fig. [2] In all panels the images are 8" by 
8", with orientation north up and east left. The left panels show 
interpolated images of the QSO at wavelengths near to the DLA 



line. Contours correspond to an image centered on the visually 
detected emission feature close to the DLA redshift. In the mid- 
dle panels in Fig. [2] the plots are reversed, such that the image 
shows the emission line object and the contours correspond to 
the QSO narrow-band image. Here, the innermost contour cor- 
responds to the seeing FWHM. To enhance the visibility of the 
candidates the QSO emission was subtracted from the data cubes 
before creating the images. This subtraction of the QSO emis- 
sion w as done using a simple approach (see IChristensen et al.l 
2006). A scale factor was determined for each spaxel by dividing 
each spectrum by the extracted one-dimensional QSO spectrum. 
Using this scale factor, the QSO emission was subtracted, a pro- 
cess which retains objects with spectral characteristics different 
from the QSO in the data cube. 

The spectra of the candidates are shown in the right hand 
column in Fig. [2] These are created by co-adding spectra from 
between 4 and 10 spaxels. The dotted line corresponds to the 
lcr noise level determined from a statistical analysis of the pixel 
values in the data cube, while the lower sub-panels show the 
background noise spectra in the data cubes, obtained from 4-10 
background spaxels. 

Properties of the candidate objects corresponding to those 
with spectra in Fig. [2] are listed in Table [3] Offsets in RA, DEC 
from the QSO and the corresponding projected distance at the 
DLA redshift are listed in columns 2, 3, and 4. Emission lines 
were fit using ngaussfit in IRAF, redshifts listed in column 5 
are derived. Fluxes in column 6 are derived from the Gaussian 
fits, and errors in the peak intensity, line width, and contin- 
uum placement are propagated to calculate the uncertainties. The 
fluxes have not been corrected for the Galactic extinction. The 
flux measurements and the associated errors indicate that most 
of the candidates are detected with a signal-to-noise ratio < 4cr. 
Column 7 gives the velocity difference between the systemic 
DLA redshift and the candidate Lya emission lines. We integrate 
the signal-to-noise estimate over the emission line (Column 8), 
S/N,„, = //( V^V x a), where / is the line flux, N is the number 
of pixels the emission line covers and <x is the noise in adjacent 
wavelength intervals. Column 9 gives the observed emission line 
FWHM after correcting for the instrumental resolution. Finally 
column 10 gives the significance classes of the candidate detec- 
tion, which is explained in Section IBTTl 

Columns 3 and 4 in Table|4]lis t the values of Galactic redden- 
ing towards each QSO (Schlegel et al. 1998), and the correction 
factors to be ap plied to the candi date fluxes for a Milky Way 
extinction curve (Fitzpatric2[T999l). 

5. Results 

This section describes the classification of the emission line can- 
didates. We estimate the contamination from spurious detections 
and from interlopers. Notes on each observed object are pre- 
sented as well. 

5.1. Candidate significance class 

To estimate how reliable the candidate detection was, various 
tests were applied to the data cubes. The candidates were as- 
signed a significance class: 1, 2, 3, and 4 according to how many 
of the following tests were passed. 

1 . Instead of co-adding all data-cubes, two independent subsets 
of exposures were created, and the emission line candidate 
was visible in both sub-set combinations. In Sect. I5.5l these 
will be referred to as subcombinations. 
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Fig. 2. Left panels: narrow-band images of the QSOs with overlayed contours of narrow -band images centered on the Lya wave- 
lengths of the DLAs. The images are 8" square. Contour levels correspond to 2, 3, 4<x levels above the background noise. Middle 
panels: the reverse, where the contours are arbitrary apart from the central one that shows the QSO seeing FWHM. Right hand 
panels: Spectra of candidate emission line objects created from co-adding spectra from spaxels associated with the emission line 
candidates. The width of the grey bars over the emission lines correspond to the wavelength ranges of the narrow band images. The 
lines below the spectra show the background sky noise spectra determined from background spaxels. 
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Fig. 2. Plots of candidates- continued. No candidate is found for the z = 2.254 DLA towards Q145 1 + 1223. 
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Fig. 2. Plots of candidates- continued. No candidate is found for the z = 3.554 DLA towards Q1802+5616. 
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Fig. 2. Plots of candidates- continued. No candidates are found for the z = 3.811 DLA towards Q1802+5616, or the z = 3.565 
sub-DLA towards Q2155 + 1358. 
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Fig. 2. Plots of candidates- continued. No candidate is found for the z = 2.51 sub-DLA towards Q2233+131. 
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Table 3. Properties of candidate Lya emission lines. Column 2, 3, and 4 list the offsets of the candidate in RA and DEC and in pro- 
jected kpc at the Lye emission redshifts (in col. 5), respectively. Column 6 lists the integrated Lya flux in units of 1CT 17 erg cirT 2 s , 
and column 7 the velocity offset from the DLA redshift. Column 8 lists the integrated signal-to-noise ratio of the Lya emission line, 
and column 9 gives the line width of the emission lines. Fluxes have not been corrected for Galactic extinction. Colum n 10 lists the 
significance class of the detec tions as described in Sect. 15. 51 'Conf.' implies candidates that were confirmed previously dlVtoller et alj 
ll998tlDiorgovski et al.|[l996h . 



2. The emission line candidate was visible in the 'simple ex- 
tractions', i.e. where a Gaussian profile was not assumed. 

3. Emission line candidates were visible in the narrow-band im- 
ages when the QSO spectrum was subtracted from the data 
cube. 

4. Emission line objects that were directly visible in individual 
or combined data cubes, or in the stacked spectra. 



In all cases, for a candidate to be considered further it was re- 
quired to be detected above 3cr in both narrow-band images and 
the associated spectra. The significance classes of the candidates 
are listed in column 10 in Table[3] and comments for each object 
are described in Section 15.51 Since the candidates were found 
from visual inspections of the data cubes, this classification was 
done to describe the candidates in a more qualitative manner. As 
the classes involve various tests on the data sets, this classifica- 



tion goes beyond the simple statistical significance in terms of to 
what cr level the object is detected. 



5.2. Non-detections 

In the data cubes where no candidates were found, we esti- 
mated the upper limits for the emission line fluxes. Spectra from 
spaxels within one seeing element (e.g. 4 spaxels corresponding 
to a seeing of 1") were co-added to create a one-dimensional 
spectrum. Artificial emission lines with varying line fluxes were 
added to this spectrum at the DLA wavelength, and Gaussian 
profile fits to these lines were used to estimate the detection level. 
The results are listed in Table|5] The varying limits are due to the 
wavelength dependent noise in the data cubes and in particular 
the presence of residuals from nearby sky emission lines. 
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Table 4. Column 2 and 3 give values of the Galactic redden- 
ing and the corresponding correction factor to be applied to the 
emission line candidates. 



We also tested an automatic routine where the re-detection 
of the artificial objects was done with no visual intervention. A 
set of narrow-band images were created in wavelength ranges 
around the artificial line. For the detection of an emission line 
the location was constrained to be within +10 A of the input 
central wavelengths. These images were smoothed and a two- 
dimensional Gaussian profile was fit to the images. When an ob- 
ject was detected above a certain threshold, spaxels around the 
centre within the seeing FWHM were co-added. A series of tests 
showed that the recovered flux was consistent within Icr errors 
for fluxes down to / = 5 x 10~ 17 erg crrT 2 s _1 . In a typical data 
cube this was also the detection limit where 50% of the objects 
were re-identified, while the fraction of re-identified emission 
lines at this flux level from a visual inspection was larger. 

Tests on the frequency of false detections in data cubes where 
no objects were present showed that simultaneous detections 
of objects in narrow-band images and associated spectra with 
S/N > 3 occurred at a rate of less than 5% in a series of ex- 
periments. Therefore false detections cannot explain the large 
number of candidate objects. 
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Table 5. DLA and sub-DLA systems where no candidate emis- 
sion lines are found and 3cr upper limits for the line fluxes. 
Fluxes are in units of 10~ 17 erg crrT 2 s _1 . 



5.3. Experiments with artificial objects 

To investigate how the efficiency of the visual inspection de- 
pended on object properties, several experiments with artificial 
data cubes were made. Similar to artificial experiments for one- 
and two-dimensional data sets, artificial emission line objects 
were added to the data cubes. These objects were described by 
the location in RA and DEC, central wavelength, peak emission 
intensity, and the widths in RA, DEC and wavelength. For sim- 
plicity we assumed that an emission line object seen as a point 
source in the data cube could be represented by a Gaussian pro- 
file in each direction, i.e. described by a Gaussian ellipsoid in 
the data cube. 

We first tested completely simulated data cubes with statisti- 
cal noise levels corresponding to the typical noise level in the 
combined data cubes. An emission line object with a flux of 
5 x 10~ 17 erg cirT 2 s , a width of 800 km s _1 , and spatial 
FWHM of 1" was placed at a previously known wavelength. In 
the stacked spectra no objects could be seen immediately. The 
emission line was only identified after inspecting the data cube 
in the visualisation tool, and it was extracted and analysed in the 
same way as the real data. Similar tests were made by adding an 
emission line to a real data cube, where the background noise in- 
cluded the systematic noise as well as the pure Poissonian noise. 
These tests produced similar results for the faint emission lines 
with Lya flux / ~ 5 x 10~ 17 erg crrT 2 s _1 , i.e. 1) the emission 
line flux could be recovered within uncertainties, 2) even at very 
small impact parameters the object could be found 3) the recon- 
structed PSF of the emission line object was irregular as in any 
of the images in Fig. [2] 



5.4. Field Lya emitters 

We estimate here whether the detected candidates are likely to 
be field Lya emitters having no association with the DLAs. 
Observations of high redshift objects have partly focused on 
detecting Lya- emission from galax ies to determine the g lobal 
comoving star-formation rates (e.g. iHu et alJU9 98. 2004). The 
density of Lya emitters at z ~ 3 is estimated to be 15000 deg~ 2 
Az _1 with line fluxes brighter than a mean of/ = 1.5 X 10" 17 

erg cm 2 s^ 1 (IHu et alj |1998t iKudritzki et al J 120 00). From the 
luminosity function at z ~ 3 (Ivan Breukelen et aLll2005l) . the ex- 
pected number of field Lya emitters at a flux limit of 5 x 10~ 17 

erg cirT 2 s _1 is 1.7xl0~ 4 arcsec~ 2 Az ■ In our survey, the 9 data 
cubes sample a total redshift interval of Az = 2 1 .55 around z ~ 3. 
Statistically, it is expected that there are 0.2 field Lya emitters in 
the whole sample presented here. Because these very faint lines 
are difficult to locate when the approximate wavelength is not 
known in advance, we did not look for field emission objects. 
The negligible number of expected field emitters furthermore 
shows that the emission candidates, if proved to be real, are un- 
likely to be interloping field Lya emitters. They are much more 
likely to be associated with the DLA galaxies. 

5.5. Notes on individual objects 

This section explains the significance of the candidates for each 
individual QSO. 

Q0151+045A. - This is a z em ~ Zabs system atz » 1 .93. After 
flux calibration, the QS O spec t rum is 2 magnitudes brighter than 
that presented in lM0ller et al.l d!998l) . The low instrument sensi- 
tivity at 3560 A combined with a variable extinction coefficient 
at Calar Alto makes the calibration uncertain. 

Extended Lya emission was observed in a regio n of 3"x6" 
aroun d the QSO mostly to the east of the QSO (Fvnb oeit al.l 
1999). Long slit spectroscopy along the long axis revealed ve- 
locity structur es of 400 km s _1 that could be interpreted as a 
rotation curve (lM<*llerlll999l ). In the IFS data extended emission 
is detected to some degree in Fig. [2] but not with the same detail 
as in the higher spa tial resolution and larger field of view data 
in lFvnbo et al.l (1 19991) . This is the only case in the sample where 
extended emission is found, but the signal is not strong enough 
to determine the velocity structure over the extended region. The 
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spectrum shown in Fig.|2]is the total spectrum co-added from the 
whole nebula. 

Q0953+4749. - This z em = 4.457 QSO has th ree DLAs at 
Zabs = 3.407, 3.891, and 4.244 dBunker et al.ll2003l) . A candidate 
associated with the lowest redshift DLA is visible in the narrow- 
band image in Fig. [2] Independent subcombinations, the simple 
extraction, and the corresponding spectra show a faint emission 
line. This emission line coincides with a sky emission line 1.6 A 
away and could be due to sky subtraction errors, so we only as- 
sign this candidate a significance class of 2. A Lya emission 
line from the DLA galaxy has been reported (A. Bunker, pri- 
vate comm.) but its line flux is below our detection limit. For the 
second DLA system at z = 3.891 the object is present in sub- 
combinations, the simple-extracted images, and in the extracted 
spectra. This candidate is assigned significance class 3. No can- 
didate is found for the highest redshift DLA to the detection limit 
reported in Table [5] 

The locations of the candidates are compared to WFPC2 im- 
ages obtained from the HST archive, but no continuum counter- 
part could be identified. 

Q1347+112. - This z em = 2.679 QSO has a DLA at z abs = 
2.471 and another possible one at z a bs = 2.05, which needs con- 
firmation from spectroscopy at higher spectral resolution. An 
emission candidate for the z = 2.471 DLA is visible in the sub- 
combinations and the extracted one-dimensional spectra. In the 
simple extraction, the spectrum has a low signal-to-noise ratio 
and the emission feature in the spectrum is faint. We assign a 
significance class of 3 to this candidate. For the z = 2.0568 DLA 
system we detect a candidate emission line object, but note an 
increase in the background noise shortwards of 3750 A. The ob- 
ject is not seen in one of the subcombinations, nor the extracted 
spectra, and therefore the candidate is assigned a significance 
class of 1 . 

A snapshot WFPC F555W image dBahcall et al.lll992l) ob- 
tained from the HST archive has a 5<x limiting magnitude of 
24.4 mag arcsec~ 2 , but no continuum counterpart can be seen at 
the location of the candidate. 

Q 1425 +606. - This z em = 3.163 QSO has a DLA at z abs = 
2.827. Because this QSO is very bright, strong residuals within 
1" from the QSO centre are present in the narrow -band image 
where the QSO emission is subtracted. A faint object offset by 
~4" to the west is visible in the narrow-band image in Fig. [2] The 
candidate is present in subcombinations and in the constructed 
spectra and is assigned a significance class of 3. In PMAS data 
cubes, spaxels in the west region are more noisy than the average 
due to an overall lower transmission. Note that the tests suggest 
a good candidate, but the impact parameter (> 30 kpc) is large. 

Q1451+1223. - This z em = 3.246 QSO has two DLAs at 
Zabs = 2.469 and z a b s = 2.254 and a sub-DLA at z a b s = 3.171. 
For the DLA system at z = 2.469 an object appears after the QSO 
subtraction close to the centre. It is caused by residuals, since the 
spectrum has no emission lines at the expected wavelength. For 
the same DLA, a region ~4" to the north west appears in both 
narrow-band imaging, subcombinations, simple extractions, and 
the constructed spectra. We therefore assign a significance class 
of 3 to it, but note that is has a large impact parameter (39 kpc). 
For the z = 3.171 sub-DLA an object is detected to the north. 
Narrow-band images from subcombinations, and simple extrac- 
tions show the emission line candidate, but the corresponding 
spectra have emission lines with very low signals. The candidate 
is assigned a significance class of 3. No candidate is found for 
the z = 2.254 DLA. 



A deep optical broad- band image o f the fi eld surrounding 
this QSO was obtained by ISteidel et al.l d!995l). who fo und no 
obvious candidates to the absorbers. IWarren et al.l d200ll) found 
one candidate offset by 3'.'9 to the south-west of the QSO in a 
NICMOS image, but this object is outside the field of view of 
the IFS data. An HST/STIS archive image shows that the emis- 
sion line candidates lie in regions where no continuum emitting 
counterpart is found. 

Ql 759+7539. - This z em = 3.05 QSO has a DLA at z abs = 
2.625 and a sub-DLA at z a b s = 2.91. The candidate detected 
for the DLA system in Fig. |2]lies near the northern edge of the 
field of view and can be affected by flat field errors. Although it 
is bright, the candidate is not visible in both subcombinations, 
and is therefore assigned a significance class of 2. A bright area 
1'.'8 south west of the QSO appears after the QSO emission is 
subtracted but it is likely due to residuals. It has no emission lines 
at the expected wavelength and is not considered further. The 
higher redshift sub-DLA system has an emission line candidate 
which is visible only after the QSO PSF has been subtracted 
from the final cube. However, the candidate is only visible in 
one out of two subcombinations and we assign this candidate a 
low significance class of 1. 

A NICMOS snapshot image showed no bright galaxies 
near the QSO to a lim it corresponding to an L* galaxy 
(IColbert & Malkanll2002l) . 

Q1802+5616. - This z em = 4.158 QSO has four DLAs at 
z abs = 3.391, 3.554, 3.762, and 3.81 1. The candidate for the low- 
est redshift DLA system is directly visible in the reduced and 
combined data cube when looking at the stacked spectra. The 
candidate can also be identified in individual subcombinations 
and in the simple extracted spectrum. Therefore this candidate is 
assigned the class 4. In a narrow-band image at the wavelength 
of Lya- at z = 3 .7652 there is an emission region to the south (see 
Fig.|2]i and the corresponding spectrum shows an emission fea- 
ture. However, this line is coincident with a faint sky emission 
line, so this candidate is assigned the class 2. No candidates are 
found for the other two DLA systems. 

Q21 55+1 358. - This z em = 4.256 QSO has a DLA at z abs = 
3.316 and three sub-DLAs at z abs = 3. 142, 3.565, and 4.212. The 
observations only cover Lya for the three lower redshift systems. 
IFS covering the highest redshift system has revealed a possi- 
ble fa int candidate emission line object (Francis & McD onnelll 
2006). The candidate Lya emission line associated with the DLA 
system is visible in independent subcombinations and in the 
simple extraction and is therefore assigned a high value of 3. 
Because of the partial spatial overlap with the QSO, the emis- 
sion from the QSO is subtracted to give a cleaned emission line 
object and the associated spectrum shown in Fig. [2] A candidate 
is found to the south for the z = 3.142 sub-DLA system. This 
object is visible in the simple extraction, and subcombinations, 
but only one associated spectrum shows a clearly detected emis- 
sion line. We assign a significance class of 3 to this candidate. 
No candidate is found for the z = 3.565 sub-DLA. 

Q2233+131. - This z em = 3.295 QSO has two sub-DLAs at 
Zabs = 3.153 and z a b s = 2.55 1. The galaxy respo nsible for the 
z = 3.153 DLA was found bv lSteidel etafl d 19951) . and follow- 
up spectroscopy confirm ed this by the detection of Lya emission 
dDiorgovski et al.|[l996l) . Previous IFS of this object suggeste d 
that the Lya- emission was extended (Christens en et alj 12004). 
This is not confirmed by the higher spectral resolution data in- 
cluded in this paper, although there appears to be some faint 
emission to the east of the object in Fig. [2] The new data and im- 
proved data reduction which optimises the signa l-to-noise ratio, 
confirm the line flux Lya line flux reported in Diorgovski et al.l 
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d 1996b . No candidate was found for the z - 2.551 sub-DLA sys- 
tem, consistent w ith the upper limit fr om a deeper Fabry-Perot 
imaging analysis dKulkarni et alj|2006l) . 



6. Properties of candidate DLA counterparts 

We proceed with a more detailed analysis of the properties of 
the detected candidate Lya emission lines. Only those candi- 
dates assigned values 3 and 4 are included. Of the eight good 
candidates, we reject two due to their large impact parameters 
(> 30 kpc). However, since they fulfill the criteria for good can- 
didates, they could instead belong to a brighter component in a 
group. The average redshift of all the DLAs in the whole sam- 
ple is Sample = 3.13 while that of the six remaining candidates 
is Zcand = 3.23, hence we find no preference for detection of ei- 
ther lower or higher redshift candidates. We emphasise that the 
candidates emission lines have fluxes that are detected at the 3<x 
level, but with this in mind we compare their properties with 
those of confirmed Lya emission lines from DLA galaxies. 

6. 1 . Line fluxes 

Fig.[3]shows the inferred line fluxes of the candidates as a func- 
tion of redshift. The triangles denote our candidates and square 
symbols indicate already confirmed objects from the literature 
dlVfoller & Warrenl Il993[ iDiorpyski et all Il996t iMoller et alJ 
U998tlLeibundgut & Robertson! 1999UM0iler et alJl2002Ll2004h . 
This figure shows that the line fluxes for the candidates are simi- 
lar to those for the previously confirmed ones, which have deeper 
observations and detection levels of 5-1 Oct. 

Fabry-Perot imaging studies of QSOs with DLAs have man- 
aged to reach similar or lower flux limits t han our IFS survey 
dLowenfhal et alJll995t iKulkarni et ali r2006). With their detec- 
tion limit some objects should have been detected if the Lya 
fluxes of DLA galaxies are around the level we find for the 
candidates and the confirmed objects. IFS is useful to look for 
emission lines as it allows us to adjust a posteriori the cen- 
tral wavelength, whereas in Fabry-Perot images, the emission 
line could fall at the wings of the filter where the transmis- 
sion is lower. Another advantage of IFS observations is the 
knowledge of the spatial QSO PSF as a function of wave- 
length wh ich allows a modeling and subtraction of the QSO 
emission (IWisotzki et alJ l200l Isiiichez et all I2004T) . This al- 
lows detection of emission lines even when they are superim- 
posed on the QSO. Nevertheless we do not detect emission line 
candidates closer than about 1" from the QSO possibly due to 
subtraction residuals. The fact that the confirmed objects are 
found at smaller impact parameters compared to the candidates 
(Sect. |63]l could indicate a bias. 

6.2. Velocity differences 

An anticorrelation is expected between the Lya luminosity and 
the velocity diff erence between the Lya emission line and optical 
emission lines ( We atherlev et alj |2005). The resonant nature of 
Lya causes a shift of the emission line towards slightly longer 
wavelengths where the photons can escape absorption. When a 
larger fraction of the blue part of the line profile is absorbed, 
the remaining emission line of lower luminosity will be more 
shifted in velocity compared to brighter ones. This explanation 
is supported by the st udy of Lyq emission lines from Lyman 
break galaxies (LBGs) dShapley et al.ll2003h . 
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Fig. 3. Line fluxes of Lya emission objects as function of red- 
shift, where square symbols represent already confirmed objects 
and triangles the candidates. The solid and dashed lines corre- 
spond to Lya luminosities of 10 42 and 10 43 erg s _1 , respectively. 

Fig.|4]shows the velocity differences between Lya emission 
lines and the DLA redshifts for the candidates as a function of 
the Lya- luminosity. There is no evidence for a correlation for the 
candidates. We note that the only candidate that shows a negative 
velocity offset is the best candidate in the sample; the z = 3.391 
DLA towards Q1802+5616. For the candidates we find an aver- 
age velocity difference of 300+580 km s" 1 , which is similar to 
the velocity differences measured for LBGs; |Pettini et afl fcOOl) 
find 560±410 km s _1 while a larger sample has Av = 650 km s _1 
betwe en the Lyq- emissio n line and low-ionisation absorption 
lines (Sh aplev et al.ll2003l) . In the case that DLAs are associated 
with bright galaxies we would expect to see large velocity off- 
sets too. Furthermore, as the line of sight towards the emission 
line object and the QSOs differ by 10-30 kpc, a larger velocity 
offset can be expected due to differences in kinematics within 
the host and its environment. Instead, if the DLA galaxy resides 
in a group, the velocity offset will reflect the velocity disper- 
sion in the group instead of being related to the host. In support 
of this idea, it has been shown that bright Lyman break galax- 
ies at z > 2 are surroun ded by gas extending to large distances 
dAdelberger et alj|20 05). Correlation studies have revealed that 
DLA s cluster on almost the same scale as LBGs (Cook eet al.1 
12006b . indicating that a similar amount of gas is present in their 
environments. 

Like flux-limited surveys, this IFS study selects the brightest 
emission component, and it is possible that the real absorbing 
galaxy is a fainter component in a group. 

In the case that DLA galaxies are related to rotating large 
disks, it can be expected that the velocity difference increases 
with impact parameter but Fig. [5] shows no clear correlation. In 
three of the confirmed DLA galaxies optical emission lines have 
velocity differences betw een -200 and 30 km s relative to Lya 
dWeatherlev et al l 12005). Some candidates have larger offsets, 
possibly affected more strongly by resonant scattering. 

6.3. Hi extension 

The average impact parameter of 16 kpc derived for the can- 
didates is larger than that expected by numerical simulations 
which favor impact parameters of b — 3 kpc for DLA galax- 
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Fig. 4. Velocity differences between the Lya emission lines and 
DLA redshifts as a function of the Lya luminosity. An average 
error bar for the Lya emission candidates is shown in the lower 
right corner. 
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Fig. 5. Velocity differences between the Lya emission lines and 
DLA redshifts as a function of the impact parameter. 



ies, and have fewer than 25% with b > 10 kpc at a l l redshifts 
dHaehnelt et all l2000t lOkoshi & Nagashimal 120051: iHou et al.1 
120051) . Larger DLA galaxy sizes of 10-15 kpc at 2 < z < 4 
are inferred by other simulations dGardner et alj|200lt) . A pos- 
sibility for the difference between observations and simulations 
is that the simulations assume a single disk s cenario, while DLA 
galaxies could exist in groups dHou et al.l2005l) . The real absorb- 
ing galaxy could be fainter and lie closer to the QSO line of sight 
than the detected candidate galaxy. 

An anticorrelation between N(Li i) and the distance t o the 
nearest galaxy is found in simulations (Gardner et al. 2001), but 
no analysis of this effect for observed DLA galaxies has been at- 
tempted. A trend for larger column density absorbers at smaller 
impact param eters was observ ed in a sample of DLA galax- 
ies at z < 1 dRaoetalJ 120031 At lower column densities in 
the Lya forest such an a nticorrelation has been shown to ex- 
ist dChen et al.ll 1 998l 1200 lh . Observations of the galaxies giving 
rise to Mg n absorption lines showed an anticorrelation between 
the impact parameters and column densities of both Mg n and 



Fig. 6. Angular impact parameter vs. redshift. The solid line 
corresponds to the size -luminosity relation for an L* B galaxy 
dChen & Lanze tta 2003), while the dotted and dashed lines cor- 
respond to galaxies with B band luminosities Lb = 10%L* and 
Lb = 1%L* B , respectively. Squares are objects from the literature 
and triangles the candidates from this survey. 



Hi dChurchill et al.l |2~000|) . but recent observations of a larger 
sample have indicated that the correlation is not always present 
dChurchill et al.ll2005l) . 

We here investigate whether the candidates show a similar 
anticorrelation using the impact parameters for the candidates as 
a proxy for the sizes of neutral gas envelopes around proto galax- 
ies. Assuming such a correlation is necessarily a rough approx- 
imation because large morph ological differences between indi- 
vidua l systems are expected dRao et aLll2003t IChen & Lanzettal 
120031) . Specifically, the possible presence of sub-clumps is ne- 
glected. 

6.3.1 . DLA galaxy sizes and luminosities 



To ana lyse the sizes of DLA galaxies at z < 1 . IChen & Lanzettal 
(2003) describe the extension of the neutral gas cloud associated 
with DLA galaxies as 



— = [hi\ 

R* \Uj' 



(1) 



assuming that DLA galaxies follow a Holmberg relation be- 
tween galaxy sizes and luminosities. Their fit to the observed 
DLA galaxies gives R* = 30/i~' kpc, t - 0.26*2'?f, where L* B cor- 
responds to a galaxy with M* B = -19.6 . Base d on the morpholo- 
gies and impact parameters lChen et al.l d2005l) argued that dwarf 
galaxies alone cannot represent the DLA galaxy population at 
z < 1 . In the case that the DLA galaxy population evolves from 
low luminosity objects at high redshifts to higher luminosity at 
lower redshifts, this will affect the expected impact parameters. 
In Fig.|6]the impact parameters of candidates and confirmed ob- 
jects are shown as a function of their redshifts. Overlayed on this 
figure are curves for the size-luminosity relation derived for low 
redshift DLA galaxies. If DLA galaxies at high redshift follow 
the low redshift scaling relation, they comprise a mix of galaxy 
luminosities. 
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Fig. 7. Column density of neutral hydrogen as a function of 
impact parameters for the candidates and previously confirmed 
objects. The symbols are similar to the previous figures. The 
solid and dashed lines are fits to the power-law relation b/b* = 
(N/N*~f for the candidates and confirmed objects, respectively. 



Fig. 8. Column density of neutral hydrogen as a function of im- 
pact parameter for the candidates and previously confirmed ob- 
jects. The solid and dashed lines show the fit to the exponential 
relation A^(H i) = A^(H i) exp{-b/h) for the candidates and con- 
firmed objects, respectively. 



6.3.2. Powerlaw profiles 

Using similar arguments as above one could expect that there is a 
relation between the impact parameter b, and the column density 
measured for the DLA. Fig.Qshows the A^(H i) measured for the 
DLA as a function of the impact parameters in kpc. We assume 
a similar scaling relation as in Eq. ([TJ for the impact parameter 
and N(H i), i.e. 



b _(N(Hi) 
b* ~ \JV(Hi)* 



(2) 



We set log Af(Hl)* = 20.3 and the error on the measured 
impact parameter is given by the fibre size of (X'5, which corre- 
sponds to ~4 kpc. A fit of the observed impact parameters for 
the candidates gives b* = 15.9 ± 1.4 kpc, and p = -0.23 + 0.08 
as shown by the solid line, while for the confirmed objects, we 
find b* = 12.0 +3.7 kpc and p = -0.36 + 0.14 as represented by 
the dashed line. The fits to the candidates and confirmed objects 
are consistent within l<x uncertainties. 

6.3.3. Exponential profiles 

Radio observations of the 21 cm emission from Hi disks in the 
local Universe have shown that an exponential profile is a poor 
representation in the central part of disk galaxies, where the 
21 cm flux de nsity either stays constant or decreases towards 
the centre (e.g. lVerheiien & S ancisi 2001). However, at optical 
wavelengths disk galaxies are well represented by exponential 
profiles. Here we fit the impact parameter distribution by the re- 
lation 



JV(Hi) = JV(Hi) exp(-b/h) 



(3) 



where h is the scale length, A^(H i) the central column density 
of a simple exponential disk, and A^(H i) is the column density 
measured for the DLAs. The resulting fit to all the candidates 
is shown by the solid line in Fig. [8j which has logA^(Hi) = 
21.7 ±1.1 cirT 2 and h = 5.1*H kpc. This result is similar within 
the uncertainties to a fit to the confirmed objects only (logA'o = 
21.7 ± 1.1 cm" 2 and h = 4.5+^ kpc). 



Local disk galaxies have optical scale lengths ranging from 
~2 kpc to ~6 kpc, and observations of the H i profile in low sur- 
face brightness galaxi es have indicated scale lengths > 10 kpc 
(Matthew s et alj|2001l) . In contrast, higher redshift spiral galax- 
ies in the HST Ultra Deep field have sm aller optical scale lengths 
of 1.5-3 kpc dElmegreen et alj 120051) . possibly biased towards 
smaller values due to the easier detection of high surface bright- 
ness, high star formation rate regions. The question is how ex- 
tended the gaseous envelopes are around these young galaxies. 
The impact parameters of the candidates suggest that high red- 
shift DLAs reside f ar from the host g alaxy, if not in a regular 
proto-galactic disk (Wo lfe et~aT] [l986). then in a region of the 
same physical scale, possibly in merging clumps of gas sur- 
rounding the actual proto galaxy. In this picture, DLAs can be 
found far from the center of the parent galaxy. 

The two objects that were originally discarded as candidates 
due to their large impact parameters (> 30 kpc) do not follow 
the relations for either the exponential or power-law profiles. 
Including them would make the scatter around the fit substan- 
tial. 



6.4. Metallicity effects on Lya emission 

Using a space based imagin g survey and follo w up long-slit 
spectroscopic observations, Nfoll er et~ai1 (2004) found indica- 
tions for a positive metallicity-Lya luminosity relation, such 
that Lya emission was preferentially observed in higher metal- 
licity systems. They argued that this positive correlation could 
over-power the negative dust-Lya luminosity effects which 
are expected to be strong in high metallicity environments 
(Chariot &Fd3 ll993|). Studies of Lya emission from nearby 
star-forming galaxies have not revealed any correl ations between 
metallicity and Lye emission strength (Keel 2005). Differences 
in the velocity-metallicity relation between high and low redshift 
DLAs could be expl ained by higher re dshift DLAs residing in 
lower mass galaxies dLedoux et al.ll2.006h . that have fainter Lya 
emission. 

In this context we investigate the distribution of metallici- 
ties for the emission line candidates in comparison to the total 
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Fig. 9. Cumulative distribution of Si metallicities of the six good 
emission line candidates compared to the remaining part of the 
sample. The probability that the two distributions are similar is 
38% estimated from a Kolmogorov-Smirnov test. 



sample. We compare the cumulative distributions of metallici- 
ties ([Si/H]) for the DLAs that have candidate Lya detections 
with the metallicities of the remaining objects that have either 
no Lya candidates or rejected candidates. The distributions in 
Fig. [9] show the fraction of DLAs with metallicities larger than a 
given value. Table[JJhas several lower limits on [Si/H] and Fig. [9] 
treats the limits as actual detections. 

A two sided Kolmogorov Smirnov (KS) test gives a prob- 
ability of 38% that the two samples have the same underlying 
distributions. A similar analysis for [Fe/H] gives the same re- 
sult. Hence, none of the tests give clear statistical evidence for a 
difference between the two populations. Only a small number of 
DLAs are included in this survey. For the two samples with Ni 
and N2 being the number of objects in each sample respectively 
we find NiN 2 /(Ni + N 2 ) = 4. For the KS test to b e statistically 
valid a value larger than 4 is required dPress et aTlll992l) . hence 
a few more objects are needed to make the test more statistically 
significant. 
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Fig. 10. DLA metallicities as a function of the impact pa- 
rameters, where symbols shapes have similar meanings as 
in the previous figures. The dotted line shows a fit to all 
the objects excluding the limits. This line has a gradient of 
-0.024 + 0.015 dex kpc 1 . Data for the confirmed o bjects are 
either [Si/H] of [Zn/H] taken from M0ll er et al.ld2004l) . 



to see a tendency for higher metallicities for the DLA galax- 
ies detected at smaller impact parameters. DLA galaxies are 
on th e average fainter th an LBGs detected in flux limited sur- 
veys (IM0ller et al.l 12002). Combining this with a high redshift 
luminosity-metallicity relation can imply low DLA metallicities 
without involving metallicity gradients. 

Fig. [10] shows the DLA metallicities as a function of the im- 
pact parameters. The line shows a fit to all objects ignoring those 
with lower limits on [Si/H]. This gradient is -0.024±0.015, i.e. 
is consistent with zero at the 2cr level. The large scatter in the 
plot could be real and unrelated to gradients. Different DLAs 
exhibit a large range of s tar formation histories dErni et al.| [2006: 
iHerbert-Fort et al. 2006), which makes it unreasonable to expect 
a smooth relation between metallicities and impact parameters 
for a sample of DLAs. Clearly more data are needed to deter- 
mine the reality of any relation. 



6.5. Metallicity gradients 



7. Conclusions 



In local galaxies the metallicities of H 11 regions ar e shown to de- 
creas e with increasing radial distance in the disk (Zaritsky et al. 
119941) . If DLAs arise in disks lower metallicities are expected 
at larger linear separations between the QSO line of sight and 
the galaxy centers. A comparison of absorption metallicities for 
3 DLAs at z < 0.6 with abundances derived from strong emis- 
sion line diagnostics from the galaxy spectra revealed that gra- 
dients are likely present at at level -0.041 ±0.012 dex kpc 1 
dChen et alj 120051) . Uncertainties in the gradients arise due to 
an unknown correction for dust depletion and also the inclina- 
tion of the galaxy plays an important role due to projection ef- 
fects. Metallicity differences due to differential depletion within 
a singe DLA galaxy c an be strong as de monstrated by observa- 
tions of a lensed QSO dLopez et al.ll2005l) . 

Observations of nebular emission lines from seven galax- 
ies at 2.0 < z < 2.5 indicate an average solar metallicity 
(Sha plev et al.ll2004j). with evidence for the presence of metal- 
licity gradients (Forster Schreiber et alj2006h . If metallicity gra- 
dients exist for high redshift DLA galaxies, we would expect 



We have presented an integral field spectroscopic survey of 9 
high redshift QSOs, which have a total of 14 DLA systems and 
8 sub-DLA systems. We detect eight good candidates for Lya 
emission lines from DLA and sub-DLA galaxies. Two of these 
are found at impact parameters larger than 30 kpc, and are not 
likely associated directly with the absorbing galaxy, but could 
be associated with galaxy groups in which the real absorbing 
galaxy resides. All candidates are detected at a statistically sig- 
nificant level in reconstructed narrow-band images as well as 
in the co-added one-dimensional spectra. Further observations 
will be useful to independently confirm the candidates at an even 
higher signal to noise ratio. We compare the properties inferred 
from the IFS data with those for previously spectroscopically 
confirmed Lya emission lines from DLA galaxies reported in 
the literature. We find that line luminosities are similar to those 
of previously confirmed objects, that the average impact param- 
eters are larger by a factor of ~2, and that some candidates have 
larger velocity offsets between the Lya emission line and the 
systemic redshift of the DLA system. 
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We analyse the distribution of DLA column densities as a 
function of impact parameters. Assuming that the average DLA 
galaxy is similar to a disk galaxy with an exponential profile, we 
show that it has a scale length of 5 kpc. Such a scale length is 
similar to disk scale lengths found for local spiral galaxies. This 
could imply that DLAs belong to large disks even a t high red- 
shifts as originally suggested by Wol fe et al.l (fl986). However, 
it is probably too simplistic to expect that high redshift DLAs 
reside in regular disks with similar structure to large local disks. 
DLA syst ems are generally not a s sociated with lumi nous galax- 
ies (e.g. IColbert & Malkanll200l iJyfoller et al.ll2002h . The large 
impact parameters found for the candidates indicate that the dis- 
tribution of H i clouds in DLA galaxies extends significantly be- 
yond the optical siz es of fainter d warf galaxies. 

Furthermore, Wol fe & Chenl d2006l) showed that high red- 
shift DLAs do not reside in extended disks that follow the lo- 
cal Schmidt-Kennicutt law for star formation. The IFS results 
presented here suggest that the Lya- emission is generally not ex- 
tended, and that star formation takes place at a distance of several 
kpc from the DLA. Hence, we may speculate that DLAs arise in 
the outskirts of proto-galaxies, for example in clouds of neutral 
gas around LBGs. In this case one would expect a significant 
scatter in the relation between the impact parameter and column 
density of the DLA since neutral clouds could be distributed ir- 
regularly around the galaxy. Contrary to expectations, the ob- 
jects show a small scatter around the relations in Figs. [7] and [8] 
Regardless of the distribution of neutral gas in DLA galaxies we 
conclude that there is a tendency to find a lower column density 
DLA with increasing impact parameter. Extending the investiga- 
tion to include DLAs and sub-DLAs with N(Hi)> 10 196 cm -2 
this tendency emerges for both the candidates and the confirmed 
objects. 

The velocity offsets between the Lya- emission lines and the 
systemic redshifts of the DLAs are larger for half of the can- 
didates compared to the confirmed objects. This could indicate 
an origin in groups of galaxies, where the DLA resides in a less 
luminous component than the galaxy detected in Lya. To deter- 
mine whether resonant scattering affects the candidate Lya lines 
more strongly and gives rise to larger velocity offsets than for 
the confirmed objects, observations of the corresponding opti- 
cal e mission lines which ar e shifted to the near-IR are required 
(e.g. lWeatherlev et al.ll2.005b . Optical emission lines furthermore 
have the advantage of being less affected by dust absorption 
and therefore are better for estimating the star-formation rates. 
Alternatively non-resonance UV lines such as C iv could be stud- 
ied. 

This survey was carried out with IFS on a 4-m class tele- 
scope, and the signals were generally near the detection limit. To 
verify this IFS method, it is necessary to get independent, higher 
signal-to-noise ratio spectra with a larger aperture telescope to 
confirm the existence of the emission lines. 
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